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ABSTRACT 

Details of how the primordial plasma recombined and how the universe later reionized 
are currently somewhat uncertain. This uncertainty can restrict the accuracy of cos- 
mological parameter measurements from the Cosmic Microwave Background (CMB). 
More positively, future CMB data can be used to constrain the ionization history using 
observations. We first discuss how current uncertainties in the recombination history 
impact parameter constraints, and show how suitable parameterizations can be used 
to obtain unbiased parameter estimates from future data. Some parameters can be 
constrained robustly, however there is clear motivation to model recombination more 
accurately with quantified errors. We then discuss constraints on the ionization frac- 
tion binned in redshift during reionization. Perfect CMB polarization data could in 
principle distinguish different histories that have the same optical depth. We discuss 
how well the Planck satellite may be able to constrain the ionization history, and 
show the currently very weak constraints from WMAP three-year data. 

Key words: cosmology:observations - cosmology:theory - cosmic microwave back- 
ground - reionization 



1 INTRODUCTION 

We are entering the era of precision cosmology, with future 
high precision cosmic microwave background (CMB) data in 
the offing. If the physics governing the evolution of the pho- 
ton distribution can be modelled reliably, this data offers an 
almost unique opportunity to measure a large number of cos- 
mological parameters accurately and distinguish models of 
the early universe. It is widely recognized that small second 
(and higher) order effects can effect the CMB power spectra 
at the several-percent level: these include the kinetic and 
therm al Sunyaev- ZeI'dov ich (SZ) effects (iHu fc Dodelsonl 
l2002t and CMB lensing jLewis fc Chaliinorll2006ft . As well 
as the cosmological parameters of most interest, there are 
however also other uncertain parameters governing the back- 
ground evolution of the universe that can be of comparable 
importance. 

Here we focus on the ionization history, parameterized 
as the ionization fraction as a function of reds hift. The broad 
detai l s of recombinat ion are well understood iPeebleslll968t 
ll993l:lHuet"ai]|l995h . Direct recombination to the ground 
state is ineffective as it releases a high energy photon that 
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immediately ionizes another atom: it can only be important 
if the photon is sufficiently cosmologically redshifted before 
encountering another atom. The dominant mechanism is in 
fact capture to an excited state followed by a two-photon 
transition to the ground state. However there are many 
excited states in both hydrogen and helium, giving many 
possible recombination channels. Furthermore level popula- 
tions may be out of equilibrium, requiring a full multi-level 
atom e volution to calcula te the recombination rate accu- 
rately jSeager et al.ll2000r) . This can only be done to the 
extent that the different transition rates are known (or can 
be calculated), in particular excited state two-photon tran- 
sition ra tes are not known very well and are potentially im- 
portant foubrovich fc Grachevl2 005h The two-photon rates 
can also differ significantly from their empty-space value 
due to induced decay by CMB photons dChluba fc Sunvaevl 
2006). 



Most current CMB power spectrum cal culations use 
the e ffective model used in the code recfast iSeaeer et alJ 
1999). However the ionization fraction from this code 
differs by several percent from a more recent calcula- 
tion in iDubrovich fc Grachev! J2005T) that includes addi- 
tional transitions. Level splitting and other neglected ef- 
fects can also have percent-level effects on the CMB 
power spectra JChluba fc Sunvaevll20fj5 iLeung et alJl2004 
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Rubi no-Martin et al Furthermore there is no full cal- 

culation with any quantification of the error. When high 
resolution CMB data are available we will either need to be 
able to calculate the recombination history sufficiently accu- 
rately that errors can be neglected, or we will need to include 
uncertainties in any analysis. Not surprisingly assuming an 
incorrect history can give biased constraint s on the cosmo - 
logical parameters inferred from CMB data et al 
Our aim in this paper is to see whether a crude parameter- 
ization of uncertainties in recombination can be used with 
future data to give reliable parameter constraints when con- 
siderable uncertainties remain in the details of recombina- 
tion. The detailed task of improving recombination models 
and parameterizing residual uncertainties is an important 
challenge for the future. 

Once recombination has proceeded to reduce the resid- 
ual ionized density to a low level the details are no longer 
important for the CMB: the electron density is so low that 
few CMB photons are scattered. However at some point first 
collapsed objects will form, and at some later time high- 
energy photons emitted by (for example) quasars or stars 
of various populations can cause the univer se to reionize 
jLoeb fc Barkanall200ll: iGnedin fc Far] 12006ft . Current ob- 
servations of quasar absorption spectra only put a lower 
bound on reionization redshift, giving evidenc e for the first 
neutral hydrogen on our light cone at z fa 6 feecker et alJ 
120011: iFan e t al. 2006). Exactly how the ionization fraction 
evolved between the low level remaining after recombination 
and z w 6 is unknown in any detail. 

The main CMB constraint on reionization comes from 
the large scale polarization signal generated by scattering of 
the CMB quadrupole during reionization. This gives a char- 
acteristic bump in the large scale polarization power spectra 
on scales larger than the horizon size at reionization, and the 
exact details of the shape of the bump c an in principle be 
used to constrain the ion ization history dKaplinghat et all 
120031: iHu fc Holder! l2003l . Scales smaller than the horizon 
size are uniformly damped, leading to a suppression of the 
acoustic peaks of e~ 2r where r is the optical depth to 
reionization; the small scales cannot therefore be used to 
constrain the details of the history, only the total optical 
depth. Beyond linear theory there are also tiny secondary 
aniso tropies on small scales that we do not discuss further 
here (IWelledll999tlHull2000l) . 

The current constraint on the optical depth from the 
WM AP three-year pola r ization observations is r = 0.09 ± 
0.03 JPaee et all l2006t ISpergel et all l2006h . significantly 
smaller than r ~ 0.17 favoured by the one- year data 
(which had unsubtracted foreground contamination). The 
new value may be more consistent with simple reionization 
scenarios, though the lower fluctuation amplitude measured 
by the three- y ear data makes this i n terpre t ation unclear; for 
discussion see lAlvarez et alJ |2006); Popa (2006). 

Since there is currently no convincing model of the 
reionization history we attem pt to constrain it in a rela- 
tively model independent way. IHu fc Holdel d2003l) suggest 
a binned fit and performed a principal component analysis 
with Fisher matrix techniques. They showed that the first 
two to three eigenmodes should be fairly well constrained 
with future observations, but other details are effectively 
unconstrained by the CMB alone. In the second part of this 
paper we investigate how the reionization history binned in 



redshift can be constrained w ith current, perfect, and simu- 
lated Planck l|Planckll2006h data. The general constraints 
we obtain are rather weak, however distinct models of reion- 
ization can be distinguished. Unless there is a convincing 
physical model for how reionization proceeds, flexible mod- 
elling of the reionization history is required in order not to 
obtain biased constraints on the optical depth (and hence 
the amplitude of primordial fluctuations and as inferred 
from the CMB). 

Since details of reionization are only important on large 
scales, and details of recombination are only important on 
small scales, the two effects are virtually independent. We 
start by describing our parameter estimation and forecasting 
methodology, then move on to consider recombination and 
reionization separately. 



2 PARAMETER ESTIMATION AND 
FORECASTING 

We make use of standard Markov Chain Monte Carlo 
(MCMC) methods to sample from the posterior distribu- 
tion of the parameters given real or forecast data. Ex- 
tra parameters gover ning the ionization history are added 
to th e code camb j Lewis et alJ I2000T) (based on cmb- 
fast iSeliak fc Zaldarriagalll996l) ) for computing the CMB 
anisotropies. This is then used w ith a modified version of 
CosmoMC jLewis fc Bridleir2 002) to generate samples from 
the posterior distribution. 

For idealized forecasting work we assume that the 
temperature T and polarization E fields (including noise) 
are statistically isotropic and Gaussian, and take the B- 
polarization signal to be negligible. We shall assume that 
foreground sources of polarization can be accurately sub- 
tracted using multi-frequency observations. The covariance 
over realizations is given by 




where Ni is an assumed isotropic noise contribution, and we 
take the noise on E and T to be uncorrelated. 

For a given sky realization one can construct estimators 
of the Ci given by 

^ Y = 21 + 1 ^ X* m Yi m (2) 

m 

so that (C? Y ) = Cf Y + Ni XY . The likelihood for the matrix 
of the estimators C; given a theoretical matrix Ci , assuming 
Tim and Ei m are Gaussian, is then 

-21ogP(C i |C ; ) = (2Z + l){Tr[c ; Cr 1 ] + log|C,|}. (3) 

We calculate the expected log likelihood for each set of pa- 
rameters 6. For some fiducial model with parameters 80 this 
is given by (log P{6\ data)) where the average is over data 
realizations that could come from the Go model. Hence the 
distribution we sample from is the exponential of 

(logP(0|0 o )) =-i(Tr[C(e )C- 1 (e)] +log|C(0)|), 

(4) 
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where the covariance matrix now includes all the different 
I, m modes (we use 2 ^ I ^ 2000). The mean log likeli- 
hood peaks at the true model — Go, and the shape of the 
likelihood encapsulates any important degeneracies in the 
data. To this extent our method is superior to Fisher-based 
estimates that can be misleading if the posterior is signif- 
icantly non-Gaussian. Posteriors obtained using the mean 
log likelihood method are approximately the same width as 
those obtained in most actual realizations of C close to the 
fiducial model. Our method of using MCMC has the benefit 
of being immediately applicable to real data, allowing us to 
test much of the parameter estimation pipeline for consis- 
tency by forecasting. Note that even though the mean log 
likelihood peaks at the true model, if the posterior is non- 
Gaussian marginalized constraints on individual parameters 
need not peak at the true model values. 

For our fiducial models we take the best fit six 
parameter WMAP three-year concordance ACDM model 
iSpergel et alj200a) . The parameters we vary are the baryon 
density Qbh 2 , dark matter density Q c h 2 , approximate ratio 
of the sound horizon to the angular diameter distance at last 
scattering 1006 1 (from which we derive the Hubble param- 
eter Ho — 100ft km s~ 1 Mpc~ 1 ), constant scalar adiabatic 
spectral index n s and scalar adiabatic amplitude A s at k = 
0.05Mpc _1 parameterized with a flat prior on log(10 10 ^4 s ). 
When considering recombination we assume sharp reioniza- 
tion with optical depth r (fiducial value r = 0.091). We 
neglect the neutrino masses, and assume negligible tensor 
and non-adiabatic modes. Fiducial values for other model 
parameters are Q b h 2 = 0.0223, Q. c h 2 = 0.104, H = 73, 
n s = 0.955, log(10 10 A s ) = 3.02. We assume a fixed helium 
fraction of 0.24; marginalizing over uncertainties in this pa- 
rameter would be trivial, with the effect depending on what 
(if any) external constraint is assumed. 

For the Planck satellite we take a toy model with 
isotropic noise N? T = iVf B /4 = 2 x 10~ 4 ^A" 2 on large 
scale s, with an effe ctive Gaussian beam width of 7 arcmin- 
utes dPlanckl E)06l. Since our purpose here is to isolate the 
effects of reionization we shall not include complicating sec- 
ondary signals such as SZ or CMB lensing, though these 
must of course be accounted for when analysing real future 
data. The lensing effect can be included easily enough and, 
once included consistently, has little effect on the recovered 
parameters at Planck sensitivity jLewij |2005'1 . Thermal 
SZ can in principle be removed using the freq uency infor- 
matio n. Kinetic SZ is more difficult to model iZahn et alJ 
2005), though expected to be a small signal at I < 2000; fu- 
ture work is required to model it sufficiently accurately for 
reliable parameter estimation. In principle it might be nec- 
essary to model the SZ signal as a function of reionization 
parameters since the kinetic SZ signal comes from the in- 
homogeneous reionization epoch. Indeed ultimately the SZ 
signal may provide a useful constraint on the ionization his- 
tory, though here we focus on what can be learnt using only 
the linear polarization signal. 



3 RECOMBINATION 

3.1 Uncertainties in the standard recombination 

Before discussing the effect of uncertainties in the recombi- 
nation history, we will briefly review the the standard re- 
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Figure 1. The ionization fraction as a function of redshift dur- 
ing recombination. The solid line is computed using RECFAST 
iSeager et alj|200d) . the dashed line includes a model of addi- 
tional transitions from lDubrovich fc Grache\J i2005ft . The bottom 
panel shows the several-percent fractional difference. 



combination dynamics as implemented in the widely used 
RECFAST code. For full accuracy one has to do a radiative 
transfer calculation accounting for many different levels of 
the hydrogen and helium atoms. However, for simplicity, we 
follow the notation of RECFAST and describe the recombi- 
nation process for an effective 3-level atom, i.e. a two-level 
atom plus continuum. Detailed balance equations for the 
proton fraction x p — n p /nn a nd singly ionized he lium frac- 
tion shcii = riHeiiAiH l ea d to iSeager et alJll999r> : 

dx„ 



dz 



— — = I c IpHHQH — /fe(l 



-hv H2s /kT M 



(5) 



(l + KHA H n H (l-ttp)) 

X H(z){l + z)(l+%(AH+fe)n H (l-ip)) ' 



and 

dz 



H(z)(l + z) 



X [XHellXerinanel — /feel (/He — £HeIl)e 



1 + Knel AHeflH (/h C — KHeIl)e 



-hv ps /kT m - 



(1 + #Hel(A H e + /feelW/He " SHell^-" W^m) ' 

(6) 

where x £ = n e /rm is the ionization fraction, rm is 
the total density of ionized and neutral hydrogen. The 
relation between the recombination coefficients a and 
the photoionization coefficients j3 is given by j3 — 
a exp(— hv2s/kTm)- The frequencies v are 

from the characteristic wavelength of the atomic transitions 
under consideration as indicated by their indices. Note that 
the hydrogen recombination rate «h in RECFAST includes 
a fudge parameter F to effectively describe the multilevel 
atom. /He is the He/H number ratio. The terms in brack- 
ets of expression (HJ and (|HJ are from the detailed balance 
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between recombination and photoionization, while the mul- 
tiplied terms take into account redshifting of the Lyman-a 
photons for the hydrogen atom and for Hel the 2 1 p~l 1 s pho- 
tons via the Tf-factors. Further the rates of the two photon 
decays are given by Ah and Ah b respectively. The recom- 
bination and photoionization rates depend on the tempera- 
ture of the baryons and photons, where the evolution of the 
baryon temperature Tm is given by 



dT M 
dz 



3H(z)(l + z)m c c 1 + / Hc + x, 



(Tm-Tr)- 



2Tm 



(i+*r 

(7) 

where Tr is the radiation temperature, cir is the radia- 
tion constant, and the Thomson scattering cross section ot- 
For further details and the exact values of the constants 
see lSeaeer etafl <2000Ul999l) . 

We can now discuss two slightly different calculations 
of the recombination history as shown in Fig. Q One cal- 
culation uses recfast, the other is altered at the several- 
percent level by allowing for a dditional transitions. This is 
achieved using the model of iDubrovich fc Grachev! (120051) 
by modifying the two-photon Einstein coefficients A to allow 
transitions from upper levels and making them temperature 
dependent, i.e. 



A^ A + 



(2<l) h("iC- 



V2a,c)/kT M 



(8) 



where the sum runs over the upper levels, v are the particu- 
lar transition frequencies and A^ q ^ the upper level Einstein 
coefficients with Qi the statistical weights of the states. The 
K factors also have a l ess im portant modification following 
IDubrovich fc Grachev! i2005l) . This particular modification 
may not be at all accurate, but we can take the difference 
between these two results to indicate the level of current the- 
oretical uncertainty in the recombination history. Account- 
ing for more transitions generally increases the recombina- 
tion paths and hence speeds up recombination. One of the 
main effects is to slightly shift the redshift of maximum vis- 
ibility, and hence the angular scale of the CMB acoustic 
peaks. However the exact shapes of the CMB power spectra 
are sensitive to the full shape of the recombination curve: 
the small scale peak suppression and polarization are quite 
sensitive to the width of recombination. For percent-level 
accuracy in the power spectra, the ionization fraction needs 
to be known to percent-level or better through the peak of 
the visibility. 

Helium recombination is more complicated than hydro- 
gen, but also less important: the significant difference be- 
tween the Hel recombination histories only has a percent- 
level effect on the observed CMB spectrum on small scales 
due to the slightly modified diffusion damping length. 

The dashed curves in Fig. [5] show the bias on pa- 
rameter constraints from Planck if an incorrect recom- 
bination history is us e d tha t differs by as much as the 
IDubrovich fc Grachev! J2005) corrections to recfast. It 
should ther efore be a matter o f so me priority to continue 
the work of ISeaeer et al.l (2000) and Dubrovich fc Grachev! 



20061) : IChluba fc Sunvaevl (12006ft : iRubino-Martin et alJ 
2006) to include all transitions that might be important, 

and to quantify the importance of poorly known transition 

rates on the predictions. 

In this paper we investigate the use of a crude ad hoc 
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Figure 2. Forecast Planck constraints on recombination and 
cosmological parameters. There are two fiducial models, one using 
RECFAST (thin lines) and one using REC FAST with extra transitions 
following IDubrovich fc Grachevl 120051) (thick lines). Each simu- 
lated data set is analysed in two ways: dashed curves show the 
constraints using RECFAST and no additional recombination pa- 
rameters; solid lines show the constraints allowing for four extra 
effective recombination parameters, as described in the text. The 
dotted lines show that using a fixed wrong recombination model 
can give biased parameter constraints; allowing for extra param- 
eters (solid lines) broadens the errors bars but gives consistent 
results for both fiducial models. 



parameterization that can be used to obtain correct cosmo- 
logical parameter constraints given the current level of un- 
certainties. We regard this as a proof of principle rather than 
a recommendation for practice. It may be a useful guide to 
where most effort needs to be concentrated to obtain more 
accurate results in future. 

We base our calculations on the recfast code. We 
evolve these same equations, with modifications incorporat- 
ing additional parameters. In recfast the effect of out of 
equilibrium distributions lead to a faster rate of hydrogen 
recombination than in a simple analysis. This is accounted 
for by multiplying the recombination coefficient by a 'fudge 
factor' F, fixed to 1.14 in recfast to match their full multi- 
level atom results. Other authors get slightly different val- 
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ues jRubino-Martin et alJl20o3) . We take F to be a free pa- 
rameter that effectively governs the speed of the end of re- 
combination. 

At earlier times recombination is dominated by 
the two-photon transitions, accounted for in recfast 
only from the lowest excit ed state. The calculation of 
iDubrovich fc Grachev! l)2005h includes additional transitions 
from higher s tates, and there are also c orrections from in- 
duced decay IjChluba fc Sunvaevl 1200(f) . We parameterize 
the uncertainties governing the two-photon transitions by 
using an effective two-photon rate Ah parameterized by 



Ah = A 



(2) 



l + R H (x e = 0.75) 



R H (x e = 0.25) 
R H (x e = 0.75) 



2(0.75-21(0 



(9) 

Here A^ is the (fixed) empty-space 2s-ls two-photon rate, 
and the two parameters Rh(x b = 0.75) and Rh(x £ = 0.25) 
measure an additional early and late contribution to the 
two-photon rate from extra transitions and induced decay. 
This model is not well motivated physically, but has the 
benefit of being fairly well constrained from the data. The 
exponential form allows it to capture exponential behaviour 
from Boltzmann-factors governing populations of excited 
states. Since the CMB is only weakly sensitive to the He- 
lium recombination, we use the single parameter 7?h c so that 
the effective two-photon Helium transition rate is given by 
Aho = A^ifec. We do not change anything else. Although 
rather ad hoc, we expect that these parameters should be 
able to describe most changes to the recombination rate, and 
have a physical interpretation as measuring an effective rate 
at four different levels of ionization. Alternative parameter- 
izations are discussed bv lHannestad fc Scherreil i200ll) . 

Fig. |5] shows the parameter constraints from our 
Planck mean log likelihood forecasting function, using 
two different fiducial models of recombination. Including 
marginalization over the four additional recombination pa- 
rameters allows the cosmological parameters to be extracted 
correctly in both cases, though with increased error bars. 
The posterior of the recombination parameters clearly shows 
that Rh{x £ = 0.75) and F can be constrained away from 
zero, though the other two parameters are poorly con- 
strained. With good polarization data the parameter F is 
in fact very benign: varying only F recovers cosmological 
parameters with essentially the same error bars as fixing 
it to 1.14 (the recfast value). However varying the rate 
of the two-photon hydrogen transitions shifts the redshift of 
maximum visibility, hence the large increase in the posterior 
range of 6 that governs the positions of the acoustic peaks. 
Since our parameterization only allows for speeding up of 
transitions relative to recfast, the posterior distributions 
are non-Gaussian with a tail to larger 9 corresponding to 
faster recombination, and a one-sided increase to the mat- 
ter density posterior. 

We have only considered basic 6-parameter cosmologi- 
cal models here, and more general models may be more de- 
generate with the recombination parameters, though a run- 
ning spectral index can be constrained well even with the 
extra parameters. 

Since details of recombination only significantly affect 
the small scale power spectrum, the impact of standard re- 
combination uncertainties is rather mild when considering 
only the current WMAP three-year data: parameter con- 
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Figure 3. Ionization history for standard RECFAST (fudge pa- 
rameter F = 1.14, solid), RECFAST with F = 1 (dash-dotted), and 
RECFAST with F = 1.14 and an input of 3 X 10 -24 eVs -1 per pro- 
ton to^av_^ojn_JiojnogeneOTis_darkjiiatter annihilation (dashed; 
seeJPadmanabhan & Finkbcincr 1200^), equivalent to F26 = 0.06 
of lZhang et alj 1200&I ). 




Figure 4. WMAP (solid line) and simulated PLANCK (dashed 
line) constraints on the parameter as a function of redshift 
2. Shading corresponds to the marginalized probability for the 
WMAP constraint and contours are at 68% and 95%. 



straints are affected only by a fraction of the error bar, com- 
parable to the effects of different priors, likelihood modelling 
approximations, and secondary signals. For this reason we 
have not presented constraints on the recombination param- 
eter space for WMAP; recfast is sufficiently accurate to 
obtain reliable parameters at the moment. 

3.2 Non-standard models 

Our ad hoc parameterization may also pick up (and partly 
correct for) many non-standard processes, for example an- 
nihilating dark matter could inject energy during and after 
recombination. Low values for the free fudge parameter F 
would be likely to fit these models better, corresponding 
to the increased ionization f raction at the end of recom- 
bination in these scenarios l|Che n fc Kamionkowskil 120041 ; 
iPadmanabhan fc Finkbeinedl200St IZhang et al.ll2006l) . Fig- 
ure shows an ionization histor y with dark matter anni- 
hilation from IZhang et alJ i2006l) , compared to a recfast 
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history with lower value of the fudge parameter F. The 
shape at the tail end of recombination is somewhat simi- 
lar, even though the (very small) ionization fraction at late 
times is significantly different. The resulting C; agree to be- 
low a percent except for few-percent change to the large- 
scale polarization signal. As we have seen F can be con- 
strained well, largely as a result of Planck's good polar- 
ization sensitivity, so models with significant annihilation 
should also be clearly distinguishable. If we treat the an- 
nihilation rate as a free parameter, but fix recombination 
otherwise to recfast, our forecast for Planck suggests ho- 
mogeneous rates today of > 3 x 10 _24 eVs _1 per proton will 
be detectable at 95%-confidence , consis tent with the forecast 
of IPadmanabhan fc FinkbeineJ (12005th This is also consis- 
tent with the constraint on the fudge parameter F shown in 
Figure Hwhen comparing the annihilation history with the 
similar F = 1 model shown in Figure [3] 

If the annihilation rate is much higher than that shown 
in Fig. [3] the fudge parameter F will no longer be able to 
mimic the decay well. This is because higher annihilation 
rates give a significant contribution to the optical depth 
from the lower redshift ionization fraction, and also change 
the I < 100 polarization signal more radically. Such large 
annihilation rates should however be clearly distinguishable 
using an appropriate model. 

Physically motivated dark matter candidates are gen- 
erally expected have a small effec t on the CMB powe r spec- 
tra due to decay or annihilation (IMapelli et al]l200rj) . Con- 
straints on dark matter annihilation from rec ent data via 
the eff e ct on the ioniz a tion h istory are given in lZhang et alJ 
j2006tk IMapelli et all l l2006() . so we do not consider non- 
standard models any further here. 

Ad hoc models for the injection of Lyman-a pho- 
tons during recombina tion could also delay the time of 
maximum visibility jPeebles et al.l200ltlDoroshkevich et alJ 
2003; iBean et alj l2003l. so more general parameterizations 
allowing for slower recombination would be needed to pick 
up and account for more general models. Models which in- 
crease the ionization fraction significantly at quite low red- 
shift c an give rise to significant optical depth ( [Peebles et alJ 
1200 ll: iDoroshkevich et alJ 12001 IBean et alJ 12003ft . The"se 
models are likely to be ultimately quite well constrained by 
the large scale polarization signal as part of the recombina- 
tion component. 

An alternative way of modifying the precise details of 
recombination is to increase the temperature of the mat- 
ter (baryons), Tm, by the injection of energy, assuming the 
production of Lyman-a photons to be ineffic ient. Such a 
possib ility has been considered previously in IWeller et alJ 
( 1999) and achieves modifications to the ionization history 
by shifting the balance between recombination and cooling. 
A simple phenomenological model would be to consider the 
input of energy with a Gaussian profile centred around z 
with width p. This can be achieved by inclusion of an extra 
term in recfast (Eq. Q) 



dT M 
dz 



Th 2 
— — \ ~ ex P 

p \ TV 



2{z-z) 2 



(10) 



which, in the absence of any cooling processes and assum- 
ing that z » p, would lead to an increase in the matter 
temperature ATm = Th. This would require the input of 
a total energy of 8.6eV(Th/10 5 K) per baryon over the pe- 
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Figure 5. CMB anisotropies and ionization histories for the fidu- 
cial models used in this analysis: sharp (solid), dashed (double 
reionization) and dragged (dash-dotted). The total optical depth 
for all three models is r = 0.1. The top left panel is the ionization 
fraction x e , top right is the visibility (with respect to z). Bot- 
tom left are the temperature-polarization cross-correlation power 
spectra, bottom right the i?-polarization power spectra. Error 
bars on the bottom plots show the noise plus cosmic variance from 
our simple model of Planck in the sharp reionization model. 



riod of heating. In any given scenario cooling processes will 
suppress the effect substantially, but it is still possible to 
increase the temperature of the IGM enough for it to have 
an observable effect on the CMB anisotropies. 

We have investigated the constraint imposed by the cur- 
rent WMAP data on Th as a function of z, marginalizing over 
0.05 < p/z < 0.3 and the results are presented in Fig. 0] 
We see that for late [z < 30) and early times (z > 1100) 
large amounts of energy (Th > 10 7 K) can be input without 
changing the CMB anisotropies; at low redshift heating ef- 
fects are swamped by reionization. Large heat inputs may 
be constrained in other ways, for example, by spectral dis- 
tortions in the CMB, however any reasonable heat input 
is not co nstrained by obser vational limits on spectral dis- 
tortions l|Weller et al.lll99fl|h There are tight constraints of 
T h < 10 4 K for 30 < z < 100, T h < 10 5 K for 100 < z < 500 
and T h < 10 6 K for z ~ 1000 (all limits are 2a). 

By considering a fiducial model with Th = 0, we have 
also used simulated Planck likelihood functions to show 
that more stringent constraints will be possible in the near 
future in the range 200 < z < 1500, as show by the dashed 
lines in Fig. 0] 
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Figure 6. Planck optical depth and amplitude constraints from 
the sharp model analysed using the sharp model (thin solid), the 
incorrect result from analysing a double rcionization model using 
a sharp model (thin dashed), and the consistent result from the 
double reionization (thick dashed) and sharp (thick solid) models 
using a binned reconstruction. 



4 REIONIZATION 

4.1 Binning the ionization fraction 

In order to constrain the reionization history of the Universe 
we bin the ionization fraction in redshift bins with 



— < z < Zl + — . (11) 



Az 

2 " ~ ' ~2~ 

We further introduce a maximum redshift z max above we 
follow the standard recombination history without reioniza- 
tion and a minimum redshift z m j n below which we assume 
complete reionization. In practice we require the ionization 
history to be smooth, so we join the bins using a tanh- 
function. We neglect Helium reionization at z ~ 3 as this 
only has a very small effect on the CMB. Note that since 
a constant x e electron density falls off as n e oc (1 + z) 3 in 
matter domination, the visibility per redshift interval scales 
like (1 + z) 1 / 2 , and hence x £ bins of fixed redshift width 
contribute somewhat more to the optical depth at higher 
redshift. The choice of best bin widths and offset is not ob- 
vious without any clear theoretical priors on the expected 
reionization history. These could be taken as additional pa- 
rameters, but for simplicity we fix them here; we choose the 
first bin to start at z ~ 6 where we think the ionization frac- 
tion first was last significantly less than unity. If the binning 
is in phase with any features in the ionization history the 
reconstruction will be much cleaner. It may therefore be a 
good idea to use two separate reconstructions offset by half 
a bin for comparison. 

We consider three fiducial reionization histories that 
all have an optical depth r « 0.1. The first is a "stan- 
dard" reionization history with sharp complete reionization 
at a redshift z ~ 12. The second scenario is a dragged out 
reionization and the third a double reionizat ion scenari o 
(for possible physical mod els and discussion see lCenl i2003l) : 
iFurlanetto fe Loebl (2005)). The latter two are modelled us- 
ing a series of narrow bins in redshift adjusted to give the 
same total optical depth. Other cosmological parameters are 
set to the values for the b est fit six parameter WMAP three- 
year concordance model JSpereel et alJl200rJl . 

Figure JSj shows the reionization histories and polariza- 
tion power spectra for the three fiducial models. The temper- 



ature power spectra are virtually identical, so polarization 
information is essential. The accuracy with which the spec- 
tra can be measured is limited by cosmic variance on these 
large scales even if noise and foreground uncertainties were 
negligible. This limits the amount of information that can 
be extracted from the CMB even with perfect data. The fast 
and dragged reionization histories give rather similar power 
spectra, however the double reionization history, which has 
two distinct peaks in the visibility function, gives a clearly 
differe nt prediction (for further discussion see lHu fc Holder] 
(2003)). If fast reionization is assumed, the inferred optical 
depth would be incorrect i f in reality there is double reion- 
ization jHolder et alj r2003). as shown in Figure |S| Since the 
small scale Ci scales with the amplitude as A s e~ 2r , an incor- 
rect optical depth also means that the CMB constraint on 
the amplitude (and hence e.g. erg) would also be incorrect. 
Modelling the reionization history therefore has the joint 
aims of learning as much as we can about how the universe 
reionized, at the same time as allowing other cosmological 
parameters to be constrained reliably. 



4.2 Binning priors 

It is important to be aware of ones priors when using dif- 
ferent parameterizations. In particular, when binning some 
function, introducing many new parameters with flat am- 
plitude priors can correspond to unintended priors on other 
parameters. For example, consider the case when we some- 
how measure only the total optical depth r. When r is the 
only parameter we might get r = to ± <r, and for simplicity 
suppose the posterior is Gaussian. Now say we bin the reion- 
ization history in N bins that contribute equally to the op- 
tical depth so that r = ~^2 n T n . Then if each bin posterior is 
uncorrelated and Gaussian we have r n = tq/N ± a/vN : we 
have measured each bin to have about the same small con- 
tribution to the optical depth without having any relevant 
data! This reflects the much larger parameter space volume 
where all the bins have approximately the same contribution 
than when most of the contribution is coming from just a 
few bins. 

Looking at it the other way, consider N bins with total 
optical depth r = J2„ T n where < r n < = r n (x c = 1), 
so we treat the contributions to r as though independent. 
If we have a flat prior on x e in each bin with ^ x e ^ 1, 
we have a flat prior on the optical depth from each bin of 



P(t„) = l/r,Y> for < r n < t,Y> (and zero otherwise). The 
prior on the total r is then 

P(t) = J dr . . . J dr N 5 - Tnj ]J P(T n ) 



N 

aK ikr TT 
e 11 

n=l 

N 



27T 



dT n P(r n )e lk 



dk ikT 
2^ 6 



n 



(12) 

The exact integral gives a prior for r that is piecewise 
polynomial, peaking at f = X/_7n /2 and with variance 
a 2 = ^n 1 "™ 1 ' 2 /!^- This approximates a Gaussian for large 
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Figure 7. Rcionization fraction from WMAP 3-year data using 
seven bins with Az = 3 (thin lines) and three bins with Az = 4 
(thick lines, x e assumed zero at z > 18). Vertical error bars show 
68% confidence regions. The results are largely prior driven, but 
bin correlations are < 0.2. 



N. For bins with equal t^ 1 ' the exact result is 



{r/r m -n) N {-l) n N 



-d) 



2(N 



HT)!' 



(13) 



Assuming we are in the region where the posterior favours 
t < f, the binning prior will therefore favour higher val- 
ues of t than using a flat prior on r. For r close to zero the 
prior goes like t n ~ 1 : you are very unlikely to get all the bins 
very close to zero. A good data constraint generally has an 
exponential likelihood function, so the prior is only logarith- 
mically important compared to the log likelihood, however 
for large numbers of bins the prior can have a significant 
effect. In the absence of a strong constraint from the data, 
assuming a flat prior on a large set of reionization bins will 
therefore favour posterior constraints where each bin has a 
small contribution to the optical depth and the total r is 
higher than you would get from no binning. It is important 
to use as much prior information as possible to get relevant 
results: If you use a prior you don't believe, you shouldn't in 
general believe the posterior either. 

4.3 Results 

Figure |7| shows the constraint from cur rent WMAP data 
jHinshaw et alJ 120061 ; iPaee et al.l 120061) when three and 
seven reionization bins are used. The WMAP polarization 
results are noise dominated, so the constraint is very weak. 
As shown in the figure, the conclusions one might come to 
depend on the prior: using more bins indicates that the ion- 
ization fraction is constrained to be lower at intermediate 
redshift than when using only three, if one considers the al- 
lowed 1 — a regions. With only WMAP data the result is 
largely prior driven. Since most people's prior does not cor- 
respond to a random set of amplitudes in seven bins (e.g. 
most people probably expect at most one local maximum), 
most people should not believe the result. When better data 
is used the effect of an unbelievable prior is much less impor- 
tant, but can still pull results in a direction you don't really 



PJ 

free z max is described in lSraergel et alJ l|2006|) , though this is 
not sig nificantly c onstra ined by the data either (as expected 
jKaplinghat et al-lEooli) ). 

In order to test how well CMB observations can con- 
strain reionization in principle we assume that all cosmolog- 
ical parameters, apart from the bins in the ionization frac- 
tion, are fixed. We then use the mean log likelihood from 
noise-free full-sky data to forecast the reconstruction con- 
straints, using data from / ^ 100 where the reionization sig- 
nal is important. We use z m i n = 7, z max = 25 and Az — 3. In 
the top panel of Figure |5] we show the reconstructed ioniza- 
tion histories for the three fiducial models. From these plots 
we see that for perfect data we can distinguish the double 
reionization scenario from the other two fiducial models at 
many sigma. Note however that the error bars are corre- 
lated, and any detailed hypothesis test should take this into 
account or work directly from the chains. Models can be dis- 
tinguished even if they appear consistent in a marginalized 
error-bar plot: the regions in the full n- dimensional param- 
eter space may be different even if the projections are the 
same. A reionization binning is a nice way to see what sharp 
redshift information is available, but models that differ by 
broader features (for example the sharp and dragged mod- 
els) may be better modelled using a more targeted parame- 
terization — for example a couple of parameters governing 
the slope of x e as it goes from x e — 1 at z ~ 6 to zero at 
higher redshift. 

We next include noise expected from the Planck satel- 
lite and allow the other cosmological parameters to vary. The 
result is shown in the middle panel of Figure |S| It is hard 
to distinguish the three fiducial models at the 2-a level by 
comparing a number of marginalized reionization redshift 
bin constraints. This is largely because there is a degener- 
acy between the bins over the reionization peak: Planck 
cannot resolve the redshift of the start of reionization accu- 
rately enough to distinguish high x e followed by low x e from 
the case where two bins contribute more equally. This is ap- 
parent in the (anti-)correlation of the error bars: when x e is 
significant there is a ~ 0.3 anti-correlation between adjacent 
points. If we include a prior that x e = at 2 > 18, the con- 
straint using Az = 4 bins looks much better (bottom panel 
of Figure [SJ, though the bins are now fine-tuned to a priori 
knowledge about the expected form of the history. Bin anti- 
correlations are still large at ~ 0.5-0.7. Our conclusions are 
not changed significantly if we instead fix the cosmological 
parameters: Planck constrain the other parameters well, 
and on large scales the reconstruction is noise and cosmic 
variance limited. 

The constraint on r when seven bins are used is con- 
sistent with the fiducial value, whichever ionization history 
is used: see Figure [(J Note that our analysis is assuming 
statistical isotropy; if the large scale CMB is not isotropic 
as suggested by numerous analyses of the WMAP data, it 
may be much harder to model or learn anything about the 
late time ionization history. Alternatively it is possible that 
the signal from a well understood ionization history could 
be used to learn more about the structure of the large scale 
universe. 

An alternative to our binned parameterization would be 
to use a reduced set of well constra i ned F isher PCA com- 
ponents identified bv iHu fc Holder! i2003T) . However these 
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Figure 8. Reconstructed reionization bins from perfect data (top) and Planck (middle and bottom) for the three fiducial models. Points 
show maximum sample likelihood positions, vertical error bars are the marginalized 95% confidence intervals, and the horizontal bars 
show the width of the redshift bins (Az = 3 top two plots, Az = 4 for bottom plot). The lines show the fiducial input model ionization 
history in each case. The bottom panel shows a Planck reconstruction for three bins with a prior that x e = at z > 18. Adjacent bins 
away from zero are correlated at the 0.3-0.7 level. 



have unphysical negative ionization fractions, and the phys- 
ical interpretation is clearer in easily distinguishable models 
using our direct redshift binning. Also the PCA components 
depend on the fiducial model (or an iterative scheme), so a 
one-off binned parameter run is simpler to use in practice. 
We have investigated the extraction of PCA components 
from our posterior bin parameter chains, however the results 
we obtained appear not to be very useful and rather dif - 
ferent from the Fisher components of lHu fc Holder! {2003). 
This can be partly explained because infinitesimal deriva- 
tives about a fiducial model generally give different covari- 
ance estimates than samples from posteriors generated using 
physical priors (e.g. that ^ x e 1) accounting for the full 
non- Gaussian posterior shape. 



5 CONCLUSIONS 

Modelling the ionization history is crucial to interpret high 
precision observations of the CMB correctly. Uncertain- 
ties in the details of recombination affect the small scale 
anisotropies, and an incorrect model can lead to strongly 
biased parameter constraints. We have made a first step to- 
wards modelling recombination uncertainties using an ad 
hoc parameterization governing the recombination rate at 
four different ionization fractions. This was good enough to 
recover parameter estimates that are consistent, but at the 
cost of larger error bars. There is clear motivation for future 
work to pin down the recombination history in more de- 
tail theoretically, with quantification of errors on any poorly 
measured important parameters. In principle future work 
could do a full multi-atom calculation for each cosmological 
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model, including free parameters with error bars for all rates 
that are not known very accurately. In practice it is likely 
to be possible to devise an accurate small set of effective 
equations that can be evolved quickly, along the lines of the 
current recfast code. Using a free parameterization has the 
advantage of being able to look for any surprises, for exam- 
ple energy injection from annihilating particles, however if 
the standard model predictions are not well nailed down it 
may be hard to distinguish subtle unexpected effects from 
uncertainties in the expected model. 

The reionization impacts the large scale CMB polar- 
ization anisotropics. Our conclusions effectively agree with 
previous work, though our direct binning approach is some- 
what different. Ideal data can do quite well at distinguish- 
ing distinct models, and even Planck may be able to sep- 
arate the most clear-cut cases. We concentrate on Planck, 
also there are ground based probes which could provide use- 
ful information about reionization. In particular the Back- 
ground Imaging of Cosmic Extragalactic Polarization (BI- 
CEP) has in principle the ability to measure features in the 
first trough in the E-mode anisotropics, which at least to 
some d egree could constrain the details of the reionization 
history (iKeating et al.l2003T) . However the error bars remain 
quite large even with ideal data, and the CMB is ultimately 
not going to be the best way to learn about the reioniza- 
tion history observationally: 21cm emission is likely to be a 
much better tracer. Allowing for different models in a CMB 
analysis is however crucial to obtain unbiased constraints 
on the optical depth (and hence the underlying perturba- 
tion variance) if the optical depth is quite large. Of course 
if the optical depth is observed to be low, so we know reion- 
ization must all have happened around z ~ 6, the modelling 
for the history becomes less important as there is then little 
room more complicated reionization scenarios. Alternatively 
a consensus may arise that th e reionization history is e x- 
pected to be monotonic (e.g. see lFurlanetto fc Loebl (|2005^ . 
in which case modelling is less important and significantly 
simpler. 
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